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ABSTRACT 


The  differentiation  matrix  for  a  spline-based  wavelet  basis  will  be  constructed. 
Given  an  n-th  order  spline  basis  it  will  be  proven  that  the  differentiation  matrix  is 
accurate  of  order  2n  +  2  when  periodic  boundary  conditions  are  assumed.  This  high 
accuracy,  or  superconvergence,  is  lost  when  the  boundary  conditions  are  no  longer 
periodic.  Furthermore,  it  will  be  shown  that  spline-based  bases  generate  a  class  of 


compact  finite  difference  schemes. 
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1  Introduction 

The  use  of  wavelets  as  a  basis  set  for  the  numerical  solution  of  partial  differential 
equations  (PDE’s)  is  a  very  active  area  of  research.  Wavelets,  Daubechies  or  Spline- 
based,  provide  a  very  convenient  structure  for  dividing  information  according  to  the 
frequency  of. the  information  at  a  given  location.  This  splitting  of  the  data  by  pro¬ 
jecting  onto  a  wavelet  basis  seems  ideal  for  the  solution  of  nonlinear  PDE’s  where 
shocks  or  high  frequencies  might  exist  in  a  small  portion  of  the  domain,  whereas  in 
the  remainder  of  the  domain  the  solution  might  be  comprised  of  very  low  frequency 
information. 

In  this  paper  the  goal  is  to  understand  the  character  of  a  spline-based  wavelet 
numerical  method.  Numerical  results  for  such  a  method  can  be  found  in  [12]  and 
[6].  This  understanding  will  be  found  through  the  differentiation  matrix.  In  [4] 
the  differentiation  matrix  was  found  for  Daubechies  wavelets  with  periodic  boundary 
conditions,  and  a  very  high  degree  of  differentiation  accuracy  known  as  superconver¬ 
gence  was  proven  to  exist.  Furthermore,  it  was  seen  that  a  Daubechies-based  wavelet 
method  corresponds  to  a  finite  difference  method  with  grid  refinement  in  regions 
where  small  scale  data  is  present.  In  the  more  complicated  scenario  where  periodic 
boundary  conditions  are  no  longer  assumed,  it  was  seen  in  [5]  that  the  superconver¬ 
gence  is  lost. 

In  this  paper  the  differentiation  matrix  for  a  spline-based  wavelet  basis  with  peri¬ 
odic  and  nonperiodic  boundary  conditions  will  be  constructed.  It  will  be  proven  that 
a  spline-based  wavelet  basis,  also,  displays  superconvergence  when  periodic  boundary 
conditions  are  assumed.  Furthermore,  the  superconvergence  is  lost  when  periodicity 
is  no  longer  assumed.  In  addition,  it  will  be  seen  that  spline-based  wavelet  methods 
generate  a  class  of  compact  finite  difference  schemes. 

Recall  that  the  differentiation  matrix  is  comprised  of  three  matrices,  X>  =  C DC: 

•  The  first  matrix  C  is  the  quadrature  matrix  which  maps  from  point  values  of  a 
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function  to  the  approximate  scaling  function  coefficients  at  the  finest  scale. 

•  The  second  matrix  D  maps  from  the  scaling  function  coefficients  of  a  function 
to  the  scaling  function  coefficients  of  the  derivative  of  the  same  function.  For 
convenience,  I  will  henceforth  refer  to  this  matrix  as  the  ‘derivative  projection’ 
matrix. 

•  In  this  paper  the  third  matrix  is  the  inverse  of  the  first  matrix:  C  =  C~l . 

This  paper  is  organized  as  follows: 

1.  Introduction 

2.  Definitions  and  Framework: 

Aliasing  is  reviewed. 

3.  Multiple  Scale  Wavelet  Decompositions: 

Wavelet  methods  use  information  at  many  scales,  but  all  properties  of  wavelet- 
based  differentiation  are  dictated  by  differentiation  in  the  finest  scale  approxi¬ 
mation  subspace,  Vo. 

4.  Differentiation  in  a  Wavelet  Subspace: 

Here  we  begin  to  develop  the  relations  needed  to  build  the  differentiation  matrix. 

5.  Theoretical  Accuracy: 

In  this  section  the  theoretical  accuracy  will  be  proved.  For  a  spline-based 
wavelet  basis  using  an  n-th  order  spline  and  assuming  periodic  boundary  con¬ 
ditions,  differentiation  is  accurate  to  order  2n  +  2. 

6.  Equations  for  Finite  Dimensions: 

In  finite  dimensions  all  wavelet  operations  can  be  written  in  matrix  form.  In 
this  section,  the  matrix  equations  will  be  derived  which  will  hold  for  arbitrary 
boundary  conditions. 
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7.  Examples  with  Periodic  Boundary  Conditions: 

Explicit  differentiation  matrices  for  the  n  =  1  and  n  =  3  spline-based  wavelets 
will  be  given  for  periodic  boundary  conditions. 

8.  Examples  with  Nonperiodic  Boundary  Conditions: 

Again,  explicit  examples  are  given.  The  most  important  point  to  note  in  this 
section  is  the  loss  of  superconvergence  at  the  boundaries. 

9.  Loss  of  Superconvergence: 

The  very  important  issue  of  the  loss  of  superconvergence  when  nonperiodic 
boundary  conditions  are  imposed  is  explored. 

10.  Conclusion: 

The  role  of  aliasing  is  noted  in  determining  the  order  of  accuracy  of  T).  It  was 
shown  in  [4]  that  Daubechies- based  wavelet  methods  have  a  strong  relationship 
with  ordinary  finite  difference  methods.  In  this  paper,  a  strong  relationship 
between  spline-based  wavelet  methods  and  compact  finite  difference  schemes  has 
been  exposed.  Furthermore,  superconvergence  is  lost  when  boundary  conditions 
other  than  periodic  boundary  conditions  are  imposed. 


2  Definitions  and  Framework 


Please  note  that  the  framework,  definitions,  and  even  much  of  the  notation  in  this 
paper  come  from  a  paper  by  M.  Unser  and  A.  Aldroubi,  see  [10]. 

2.1  Fourier  Transforms 


Let  us  set  once  and  for  all  the  definitions  of  the  Fourier  transform  and  the  semi¬ 
discrete  Fourier  transform  used  in  this  paper.  The  Fourier  transform  of  f(x)  €  L2(R) 
is, 

Ho  =  r  noc^dx,  a) 

y— oo 

and  its  inverse  is, 

=  /«)*-«■<*{■  (2) 

The  semi-discrete  Fourier  transform  (SDFT)  of  fk  €  l*(Z)  is, 


and  its  inverse  is, 


/«>=  £:  ***/*, 

fc=— oo 


(3) 

(4) 


2.2  B-Splines 


All  B-splines  in  this  paper  are  central  B-splines:  centered  at  0.  The  B-spline  of 
order  n  is  denoted  by  /?n(z).  All  splines  of  order  n,  denoted  by  ■sn(x),  are  a  linear 
combination  of  B-splines  of  order  n,  see  [7]  and  [8]: 


=  (5) 

k 

and  the  space  spanned  by  splines  of  order  n  is  denoted  by  Sn.  The  B-spline  0°(x)  is 
the  box  function,  which  is  1  for  x  €  [—1/2, 1/2]  and  zero  otherwise.  All  B-splines  of 
higher  order  are  generated  from  /?°(z): 

F'(x)  =  p*p-'(z)y  (6) 
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where  V  denotes  the  convolution  operator.  Other  useful  forms  of  this  equation  are, 


(7) 

and 

= r  0"(»)/jn(*  -  »><&. 

oo 

(8) 

and  finally, 

P'+H*)  -</>*, /»%*-•)  >. 

(9) 

where  for  f,g  €  L2(R), 

<  fi9  >—  f  f(x)g(x)dx. 

«/— OO 

(10) 

We  will,  also,  need  the  samples  of  B-splines  at  the  integers: 


&"(*)  =  (11) 

In  order  to  use  B-splines  to  construct  a  multiresolution  analysis  it  is  necessary 
that  contractions  and  expansions  of  n-th  order  B-splines  also  be  contained  in  Sn. 
This  is  only  true  when  n  is  odd,  which  is,  therefore,  assumed  for  this  paper. 

2.3  Fourier  Transforms  of  B-Splines 

From  equation  (6)  we  see  that  the  Fourier  transforms  of  B-splines  are  generated 
recursively.  That  is,  the  Fourier  transform  of  a  convolution  is  the  product  of  the 
Fourier  transforms  of  the  convolved  functions,  Therefore,  the  Fourier  transform  of  a 
B-spline  is, 

/HO  =  (12) 

where  sinc(()  =  and  is  the  Fourier  transform  of  /3°(x). 

2.4  Aliasing  and  Fourier  Transforms 

The  relationship  between  the  Fourier  transform  and  the  semi-discrete  Fourier  trans¬ 
form  and  the  role  of  aliasing  is  at  the  heart  of  a  proof  that  will  be  given  later  in  this 
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paper.  Therefore,  I  will  reconstruct  the  straightforward  algebra  leading  to  aliasing 
using  notation  defined  above. 

Let  be  the  samples  of  /9*(x)  as  defined  above: 


k = /**<*)  =  ^  jT  no'-**#. 


(13) 


Break  the  integral  up  as  follows: 

1  /(Jr+i)*  »  ... 

which  is  equivalent  to, 


(14) 


«  =  i  E  r  n(  +  2xr)e-«‘if,  (15) 

2*  rs=— oo  •'”» 

where  I  have  used  the  fact  that  1  =  e~'2rrk.  If  we  compare  this  expression  with 
definition  (4)  we  see  that, 


6"(f)  =  £/»"(«  +  2xr),  (16) 

r=-oo 

which  illustrates  the  aliasing  of  frequencies  £  -f  2xr  onto  frequency 


2.5  Spline-Based  Wavelet  Bases 

As  mentioned  above,  the  subspace  of  n-th  order  splines,  5",  is  spanned  by  n-th  order 
B-splines.  The  following  theorem  will  define  the  scaling  functions  for  a  spline-based 
wavelet  basis,  see  [11]: 

Theorem:  The  set  of  functions  {^n(x  —  k) :  k  €  Z}  with 

n*)  =  E  Pkir(*  -  *)>  (17) 

fcs-00 

is  a  basis  of  5"  provided  that  the  sequence  {p}  is  an  invertible  convolution  operator 
from  1%  into  itself. 

The  sequence  {p}  is  chosen  such  that  the  scaling  function  satisfies  whatever  prop¬ 
erties  are  specified.  That  is,  one  can  specify  that  the  scaling  functions  and  wavelets 
be  orthogonal  which  places  restrictions  on  {p}.  There  are  many  ways  to  define  the 
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scaling  function,  but  in  this  paper  it  will  be  seen  that  in  the  calculation  of  the  dif¬ 
ferentiation  matrix  that  the  sequence  {p}  will  ‘divide  out’  producing  a  result  which 
holds  for  all  spline-based  wavelet  bases  produced  using  the  current  framework. 

Throughout  the  paper  I  will  use  the  notation  fik(x)  to  denote  translation  of  4>n(x) 
by  k: 

«(*)  =  Pit  ~  k).  (18) 

Define  the  dual  of  <f>n,  <f>n,  as 

*’(*)  =  £  rtP(x-k),  (19) 

k=—oo 

where  the  sequence  {r}  is  chosen  such  that 

fa  =<  to,  tk  >  •  (20) 

Later  the  relationship  between  {r}  and  {p}  will  be  made  precise.  However,  as  men¬ 
tioned  above,  {p}  and  {r}  will  have  no  affect  on  the  differentiation  matrix. 

The  wavelet  is  defined  as, 

tn(x)=  £  qk$n(x-k ),  (21) 

OO 

and  its  dual  is 

tn(x)  =  J  sfc/3n(x  -  fc),  (22) 

k=—oo 

where  $  is  the  ‘B-spline  wavelet’.  For  this  paper  this  B-spline  wavelet  is  not  needed 
and  will  not  be  dealt  with,  see  [10].  Furthermore,  as  with  the  scaling  function,  the 
sequences  {q}  and  {s}  defining  the  wavelet  and  its  dual  will  be  of  no  consequence  in 
the  calculation  of  the  differentiation  matrix. 

This  completes  the  definitions  and  framework. 
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3  Multiple  Scale  Wavelet  Decompositions 

One  of  the  underlying  strengths  of  a  wavelet  decomposition  is  that  it  is  possible  to 
easily  decompose  a  function,  or  signal,  into  its  components  at  various  scales.  To 
illustrate,  suppose  that  we  are  restricted  to  a  finite  number  of  dimensions,  say  d. 
Denote  the  finest  scale  approximation  subspace  by  Vo,  and  the  coefficients  of  the 
expansion  in  V0  by  s.  A  usual  wavelet  decomposition  would  appear  as, 


V0  =  W,  0  W2  ©  ...  ©  Wj  ®  Vj,  (23) 

where,  as  is  customary,  Wj  denotes  the  wavelet  subspace  at  scale  j.  A  function 
projected  onto  the  above  subspaces  would,  therefore,  be  represented  by  d  coefficients, 
say  c.  Under  the  assumption  of  finite  dimensionality,  there  will  be  an  invertible  matrix 
mapping  from  s  to  c: 

c  —  PS 


That  is,  if  g(x)  €  Vo  with  the  expansion, 


g(x)  =  £  s^(x), 
»=o 


and  one  projects  into  V\  and  W\  by, 


d/  2-1  d-1 

Pvx9{x)  =  EEs!<  (*)» 

j=0  i= 0 


and 


(24) 


d/  2-1  d- 1 
Pw>g(x)  =  £ 

j=0  i=0 

then  the  matrix  form  of  the  above  linear  transformations  is  denoted  P: 


(25) 


(26) 


'  Pi.‘  =<  i)  >. 

and 

Pd/2+j.i  =<  >, 


for  t  =  0,  1  and  j  =  0,  ...,d/2  —  1. 


(27) 


(28) 
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Suppose,  now,  that  we  have  found  the  matrix  D  mapping  from  s  to  i,  where  s 
represents  the  coefficients  of  the  approximation  in  V0  of  the  derivative  of  a  function. 
Of  course,  no  one  would  choose  to  work  with  wavelets  and  stop  the  decomposition 
in  the  subspace  Vq.  The  matrix  D  does,  however,  characterize  differentiation  in  any 
combination  of  wavelet  subspaces,  and  the  coefficients  c  can  be  found  from, 

c  =  PDP~xc.  (29) 

To  restate,  D  dictates  the  character  of  wavelet  differentiation.  The  primary  con¬ 
cern  in  this  paper  is  to  understand  D.  Therefore,  the  remainder  of  this  paper  will 
not  mention  wavelet  decompositions  at  multiple  scales. 
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4  Theoretical  Accuracy 

The  accuracy  of  the  spline  based  differentiation  matrix  can  be  proved  by  observing 
the  behavior  of  the  semi-discrete  Fourier  transform  of 

dq=<<j> 0,4>q>  (30) 

near  £  =  0.  It  will  be  shown  that  the  Taylor  series  of  d(()  about  £  =  0  is  of  the  form, 

d(()  =  i(  +  0(f3n+3)  +  ....  (31) 

and  this  implies  accuracy  of  order  2 n  +  2. 

4.1  The  Semi-Discrete  Fourier  Transform  of  {<f} 

As  above,  define  dq  as, 

dq=<<j>o,j>q  >,  (32) 

and  recall  from  the  previous  section  that  this  is, 

dg  =  J2pirm  /  0n(x  -  m)Pn(x  -  l  -  q)dx.  (33) 

l,m  J 

Again,  note  that  as  of  now  all  summations  and  integrals  are  from  —  oo  to  oo  allowing 
one  to  shift  summations  and  integrals  without  changing  the  limits, 

dq  =  /  P*{x)P*{x  +  m-l-  q)dx.  (34) 

Let  us  rename  this  integral  as, 

a/+,-m  =  J  /3n(x)4n(x  +  m-l- q)dx,  (35) 

which  yields  the  following  expression  for  d„ 

dq  =  £pirmai+q-m.  (36) 

l,m 

Break  apart  the  double  summation, 

dq  —  Y.Pl  y  (37) 

I  "» 
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and  shift  the  index  /, 


dg  —  ^  1  Pi— 9  y  '  rm&l—i i 


Let  the  superscript  “  '  ”  denote  time  reversal  of  a  sequence  to  get, 

d,  =  YIp'i-I  £  rmal-m-  (39) 

/  m 

We  see  that  the  right-hand  side  is  a  double  convolution.  We  can,  therefore,  use  the 
convolution  theorem  to  get, 

dtt)  =  (40) 

Recall  from  the  previous  section  that  the  following  expression  for  f(£)  was  found: 

*«)  =  (p'«)iJ”+,(or'-  hi) 

A 

Using  this  expression  for  r(()  in  the  expression  for  d(£)  we  get, 


-  pn+t((y 

A 

The  next  two  subsections  wil  derive  a(£)  and  b(£).  After  these  two  expressions 
are  found  the  proof  of  accuracy  will  begin. 

4.2  The  Semi-Discrete  Fourier  Transform  of  {a} 

Define  the  projection  of  the  derivative  of  a  B-spline  onto  a  B-spline  as, 


a(x)  =  f  /9n(y)4n(y  -  x)dy. 

«/— OO 


Recalling  that, 


m = */(*). 


and  from  the  convolution  theorem,  see  the  introduction,  we  get  the  Fourier  transform 
of  a(x), 

d(f)  =  ifW+1(0-  (45) 


The  sequence  av  is  defined  as  the  samples  of  a(x)  at  the  integers, 

*p  =  J  -  p)dx.  (46) 

The  semi-discrete  Fourier  transform  of  ap  is,  due  to  aliasing  (see  the  introduction), 

a(0  =  t  53  (£  +  2xr)stnc2n+2(f  +  2xr).  (47) 

rss— oo 

4.3  The  Semi- Discrete  Fourier  Transform  of  {6} 

Recall  that  the  samples  of  the  B-spline  /32n+1(x)  are  found  from, 

62n+l  =  j  _  k)dx.  (48) 

We  know  the  Fourier  transform  of  £2n+1(x)  from  the  introduction: 

42n+1(0  =  smc2n+2(f).  (49) 

Also,  from  the  introduction  we  know  how  to  get  the  semi-discrete  Fourier  transform 
of  bln+i  from  02n+1(f), 

62n+1(f)  =  £  W+*(f  +  2»«).  (50) 

$=— OO 

The  semi-discrete  Fourier  transforms  «(f )  and  6(f)  will  now  be  used  in  the  proof 
of  accuracy. 

4.4  Proof  of  Accuracy 

The  proof  of  superconvergence  begins  by  first  proving  two  lemmas. 

Lemma  A:  The  series, 

*i(0  =  13  smc3n+3(f  +  2xg),  (51) 

qz s— oo 

is  an  even  function  of  f  and  has  a  Taylor  series  about  f  =  0  of  the  form, 

ai(0  =  1  +  ci£2  +  +  •—  (52) 
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Proof:  Using  the  definition  of  a  sine  function  we  get, 


*»(0  = 


^  sin3n+3(£  +  2*g) 

ho  (*  +  2 *?)3"+3  • 


We  can  simplify  the  numerator  since  sin(£)  is  periodic  with  period  2x, 


~  sin3n+2(0 

,h(t+2*o2n+2' 

To  show  that  $i(£)  is  even  first  let  q  =  O  to  get, 


sin 


2n+2 


(0 


£2n+2 


=  (1  -  l/6£3  +  1/120£*  +  ...)2n+2, 


which  is  even.  Now,  sum  the  terms  for  q  =  —  1  and  q  =  1, 


9(0  = 


sin3n+3(0  sin2n+3(0 

(£  +  2x)2»+2  +  (£  —  2x)2n+2  ’ 


(53) 


(54) 


(55) 


(56) 


to  get, 

9(0  =  9(~0-  (57) 

Matching  every  term  with  positive  q  with  every  term  with  negative  q  we  see  that 
S|(4)  is  an  even  function. 

We  can  see  from  equation  (55)  that  the  first  few  terms  in  the  Taylor  series  when 
q  =  0  are  1  +  ci£3  +  c2£4  +  ....  On  the  other  hand,  when  q  ^  0  each  term  in  the 
summation  is  of  the  form, 

sin2n+2(£) 

(£  +  2x9)2"+2 

or 

(t  +  2  xV)2n+2P(0  =  sin3"+3(0.  (59) 

Using  the  series  expansions  of  these  expressions  we  get, 

(«i  +  a2(  +  a30  +  ...)P(0  =  ((-  ?/6  +  £5/120  f  ...)3n+2  (60) 


or 


(«i  +  a2£  +  a3(2  +  ...)P(0  =  O'*2  +  h.o.t.  (61) 
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This  implies  that, 

P(()  =  0(^+2)  +  h.o.t.  (62) 

That  is,  all  tenns  in  «i(0  with  q  /  0  contribute  terms  0((3n+a),  and  therefore,  all 
low  order  terms  come  from  equation  (55).  This  yields  the  desired  result.  // 


Lemma  B: 

The  series, 

»,«)  =  E  rW+a({  +  2xr), 

fs-oo 

is  an  odd  function  of  (  and  has  a  Taylor  series  about  (  =  0  of  the  form, 

MO  «  “if*1"  +  ■■■■ 


(63) 


(64) 


Proof:  First  of  all,  due  to  the  2x  periodicity  of  sin({),  sa(()  becomes, 

(£)  SL  (f  +  2rr)^- 


(65) 


Note  that  the  r  =  0  term  in  this  sum  makes  no  contribution  to  sa(£).  Now,  note  that 
each  term  in  the  sum  has  the  following  form, 

...  .  iinJ#+1(()  .  . 

/(('r)  =  r({  +  2xr)^'  (66) 

and,  consequently,  the  function  g(t  r)  defined  as  follows  is  an  odd  function  of  (: 

9it  r)  =  /tf,  r)  +  /({,  -r)  =  -(/(-*,  r)  +  f(~t  -r))  (67) 

This  implies  that  sj(^)  is  an  odd  function  of  £  since, 

5a(0  =  Y,9(M-  W 

r=l 
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Therefore,  the  Taylor  series  of  *](£)  about  {  =  0  must  contain  only  odd  powers  of  (. 
Consider,  now,  a  typical  term  of  s2(f)  for  r  ^  0, 


p(o  = 


sin 


((  +  2xr)a»+2‘ 


Rewrite  this  to  get, 


(£  +  2wr)3n+2P(()  =  sin2n+3(0. 


(69) 


(70) 


Note  that  the  binomial  on  the  left-hand  side  contains  all  powers  of  £,  and  that 
sin2n+2(f)  contains  only  even  powers  of  £  starting  with  the  term  (2w+2: 


(Co  +  cxi  +  c2e  +  =  *2n+2  +  dtp**  +  ...  (71) 

We  see  now  that  the  first  term  of  P(()  must  be  £2n+2/Co,  but  this  term  will  also 
multiply  the  cx(  term  in  the  binomial  expression  requiring  that  P(£)  have  as  a  second 
term  c*i£2n+3  where  ax  =  — ci/(c{j).  The  point  is  that  there  exist  a  nonzero  term  with 
power  2n  +  3  in  the  polynomial  P((). 

That  is,  in  general  the  series  expansion  of  the  function  /(£,  r)  about  £  =  0  contains 
all  powers  p  of  £  with  p  >  2n  +  2,  and  when  added  to  /(£,  — r)  the  even  powers  of  ( 
must  cancel  leaving  as  the  first  term  the  term  with  the  power  2n  +  3.  That  is, 

s2({)  =  of2"*3  + ...  (72) 

This  is  the  desired  result.  // 

Theorem:  The  semi-discrete  Fourier  transform  of  the  sequence  d*  has  a  series 
expansion  about  (  =  0  of  the  following  form: 

<i«)  =  «  +  o{(*+3)  + ... 


Proof:  Recall  that  the  semi- discrete  Fourier  transform  of  {d}  is, 


d(n  =  &K) 
&*»+*(()■ 


(73) 
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Now,  substitute  in  the  previously  derived  expressions  for  &({)  and  to  get. 


_  i  E£-c(l  +  2Tr)stnca*+a(£  +  2*r) 

E£-<»  +  2t9) 


f«*oo 

and  combine  similar  terms  to  get, 


flat— oo 


fa- oo 


Now,  apply  the  results  from  Lemmas  A  and  B  to  get, 

(<i({)  -  ifl(l  +  cf  + ...)  -  ai™  + .... 


This  equation  implies  the  desired  result, 


HO  =  * + ok*-*3) + ... 


This  completes  the  proof.  // 


(74) 


Multiply  both  sides  by  the  denominator  to  get, 

HO  jr  +  2*,)  =  i  £  «  +  2*r)  W«({  +  2*r),  (75) 


(d(0  -if)  Y.  +  2»l)  =  2»  E  rW+»(f  +  2*r).  (76) 


(77) 


(78) 


This  proves  that  the  sequence  {d}  differentiates  with  accuracy  2n  +  2,  and,  conse¬ 
quently,  the  differentiation  matrix  containing  the  periodic  version  of  {d}  as  its  rows 
will  differentiate  with  accuracy  2n  +  2. 

The  next  section  will  provide  explicit  examples  of  the  differentiation  matrices  for 
the  first  and  third  order  B-splines.  Also,  the  accuracy  will  be  shown  computationally. 


16 


5  Equations  for  Finite  Dimensions 

The  equations  that  will  be  derived  in  this  section  will  hold  for  any  spline-based  wavelet 
basis  defined  on  a  finite-dimensional  domain.  This  includes  the  equations  for  both 
periodic  and  nonperiodic  boundary  conditions.  Note  that  in  this  section  all  summa¬ 
tions  are  from  0  to  d  —  1,  and  consequently  the  summation  limits  will  often  not  be 
shown. 

5.1  Quadrature 

Let  g(x)  denote  the  projection  of  f(x)  €  L2(R)  in  Vo: 

g(x)  =  PvJ(x)  =  j  si<f>i(x ),  (79) 

1=0 

where 

=<  f,4t>  •  (80) 

Note  that  in  this  section  the  notation  bm(x)  will  be  used  to  denote  the  B-spline 
basis  functions.  For  periodic  boundary  conditions  one  could  have  simply  used  fi(x  — 
m),  but  for  nonperiodic  boundary  conditions  the  basis  functions  near  the  boundary 
are  truncated  B-splines.  The  notation  used  here  is,  therefore,  general  enough  for  both 
periodic  and  nonperiodic  problems.  That  is,  for  a  nonperiodic  domain  6o(x)  will  be 
only  the  right-hand  portion  with  support  width  1  of  the  B-spline  being  used.  Now, 
using  the  B-spline  expansion  of  the  scaling  function  in  equation  (79)  we  get, 

9(*)  =  ZsiZpiMx).  (81) 

l  m 

Combine  the  summations, 

9{x)  =  J2aiPi,mbm{x)y  (82) 

l,m 

and  sample  g(x)  at  the  integers  to  get, 

9(k)  =  yi-3tP>.mfrm(fc),  (83) 

l,m 
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or  amply, 


9k — yi 

Group  the  terms  as, 

(84) 

9k  — 

1  "» 

and  rename  the  inner  group  as  7  to  get, 

m 

9k  =  Wi- 
1 

Now,  transpose  7, 

(86) 

9k  = 

1 

and  utilize  matrix  and  vector  notation  to  get, 

(87) 

S=(Pbfl. 

(88) 

This  expression  is  a  relationship  between  the  samples  of  g(z)  €  Vo  and  the  exact 
scaling  function  coefficients.  In  general  we  need  a  relationship  between  the  samples 
of  f(x)  €  L2(R)  and  an  approximate  set  of  scaling  function  coefficients  <r, 

/=  {Pb)Ta.  (89) 

This  completes  the  matrix  form  of  the  quadrature  formula.  The  next  subsection  will 
derive  the  matrix  form  of  the  relationship  between  the  scaling  function  and  its  dual. 

5.2  The  Scaling  Function  and  its  Dual 

Recall  once  more  the  relationship  between  the  scaling  function  and  its  dual, 

kj  =<  >  •  (90) 

In  terms  of  the  B-spline  expansions  we  get, 

hi  =  /  h>(x)H*)dx.  (91) 
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Let  B  denote  the  matrix  form  of  the  above  integral.  That  is,  B  contains  the  projec¬ 
tions  of  the  B-splines  projected  onto  other  B-splines.  The  above  equation  becomes, 


/  *» 

(92) 

Now,  rename  the  expression  in  the  brackets  7  to  get, 

1 

(93) 

Transpose  r  to  get, 

6ij  =  52  wL- 

1 

(94) 

The  above  equation  has  the  following  matrix  form: 

/  =  PBRt 

(95) 

This  is  the  matrix  form  of  the  requirement  that  the  scaling  function  be  orthogonal 
to  its  dual  under  translation.  The  next  subsection  will  derive  the  derivative  projection 
matrix. 


5.3  The  Derivative  Projection  Matrix 

Recall  from  an  earlier  section  that  the  elements  of  the  derivative  projection  matrix 
come  from  the  following  inner  product: 

=<  >  ■  (96) 


Replace  the  scaling  functions  with  their  B-spline  expansions  to  get, 

j  t>m(x)b,{x)dx. 


m,l 


Let  the  matrix  A  contain  the  results  of  the  above  integration, 


Now,  group  the  terms  as, 


dj,i  —  n 

mj 


dj.i  — 

m  I 


(97) 


(98) 

(99) 
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and  rename  the  bracketed  expression  as  7  to  get, 


dj,i  — 


(100) 


Transpose  p  to  get, 

dj,i  =  (101) 

m 

and  the  matrix  form  of  this  equation  is, 


D  =  RAPt. 


(102) 


This  completes  the  derivation  of  the  derivative  projection  matrix.  The  following 
subsection  will  combine  the  results  from  all  the  subsections  in  this  section  to  get  the 
matrix  form  of  the  differentiation  matrix. 

5.4  The  Differentiation  Matrix 


The  differentiation  matrix  is, 


where 


Prom  above  we  found  that, 


and 


Prom  equation  (106)  we  get, 


V  =  CDC~X, 

(103) 

C  =  ( Pb)T . 

(104) 

D  =  RAPt , 

(105) 

/  =  PBRt. 

(106) 

/  =  RBtPt. 

(107) 

Solve  for  R  and  use  the  fact  the  B  is  symmetric  to  get, 


Use  this  expression  in  equation  (105)  to  get, 

D  =  (Pt)-1B~1APt.  (109) 

Now,  using  the  expression  for  the  quadrature  we  get  the  differentiation  matrix, 

V  =  bTPT(PT)-1B''APT(PT)-'(bT)-\  (110) 

or 

V  =  bTB'xA(bT)-\  (111) 

Let  me  restate  two  notational  definitions.  First,  the  elements  of  the  matrix  ‘6’  are 
bm,k  which  denotes  sampling  basis  ‘m’  at  position  ‘Jfc\  Second,  the  elements  of  the 
matrix  ‘i?’  are  /  bm(x)bk(x)dx . 

Note  that  in  the  case  of  periodic  boundary  conditions  that  all  three  of  the  ma¬ 
trices  6,  A,  and  B  are  circulant.  They,  therefore,  commute  yielding  B~x  A  as  the 
differentiation  matrix. 
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6  Examples  with  Periodic  Boundary  Conditions 

In  this  section  explicit  examples  for  the  1st  and  3rd  order  spline- based  wavelet  differ¬ 
entiation  matrices  will  be  given.  First,  the  accuracy  will  be  illustrated  explicitly  by 

^  - 

calculating  the  semi-discrete  Fourier  transforms  for  the  filters  a,  b,  and  d.  Then  the 
matrix  forms,  A ,  B,  and  D  will  be  given. 

6.1  The  Differentiation  Matrix  for  a  1-st  Order  B-Spline 

The  first  order  spline  is  the  convolution  of  two  box  functions.  It  is  defined  as  /?*(x)  = 
x  +  1  for  x  €  [—1,0],  Fix)  =  1  —  x  for  x  €  (0, 1]  and  0  otherwise. 

6.1.1  The  Filters  a,  b,  and  d 

For  the  1-st  order  B-spline  the  filters  S  and  b  are, 


o  =  {-1/2, 0,1/2), 

(112) 

and 

P  =  {1/6, 2/3, 1/6). 

(113) 

The  filter  d  is  defined  such  that, 

63  *  d  =  a. 

(114) 

Therefore, 

.3s 

II 

(115) 

where, 

a(£)  =  t  sin((), 

(116) 

and, 

63(0  =  1/3(2  +  cos(0). 

(117) 

d(£)  will  have  the  form, 

(118) 
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and  must  satisfy 

(»f  +  cC)  1/3(2  +  co«(())  =  isin«).  (119) 

Consider  the  first  few  terms  of  each  series, 

(t{  +  c£*)l/3(2  +  1  -  ?/2  +  (*/24  +  hot)  =  i((  -  e*/6  +  ^*/120  +  hot),  (120) 


or 


(t£  +  <*«)(!  -  ?/6  +  r*/72  +  hot)  =  i({  -  £3/6  +  e*/120  +  hot),  (121) 

and  finally, 


i((  ~  ^/6)  +  i(‘/72  +  cf*  +  hot  =  t(e  -  ^3/6)  +  i?/l 20  +  hot,  (122) 
with  the  conclusion  that 

9  =  5.  (123) 

This  is  exactly  the  requirement  for  4-th  order  accuracy. 


6.1.2  The  Differentiation  Matrix 


In  matrix  form,  the  filter  a  is, 


'  0 

1/2 

0 

0 

0 

0 

0 

0 

0 

-1/2 

-1/2 

0 

1/2 

0 

0 

0 

0 

0 

0 

0 

0 

-1/2 

0 

1/2 

0 

0 

0 

0 

0 

0 

0 

0 

-1/2 

0 

1/2 

0 

0 

0 

0 

0 

A  __ 

0 

0 

0 

-1/2 

0 

1/2 

0 

0 

0 

0 

A  — 

0 

0 

0 

0 

-1/2 

0 

1/2 

0 

0 

0 

0 

0 

0 

0 

0 

-1/2 

0 

1/2 

0 

0 

0 

0 

0 

0 

0 

0 

-1/2 

0 

1/2 

0 

0 

0 

0 

0 

0 

0 

0 

-1/2 

0 

1/2 

.  1/2 

0 

0 

0 

0 

0 

0 

0 

-1/2 

0 

(124) 


and  the  filter  b  is, 


'  2/3  1/60  0  0 

1/6  2/3  1/6  0  0 

0  1/6  2/3  1/6  0 

0  0  1/6  2/3  1/6 

0  0  0  1/6  2/3 

0  0  0  0  1/6 

0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

1/6  0  0  0  0 


0  0  0  0  1/6' 
0  0  0  0  0 

0  0  0  0  0 

0  0  0  0  0 

1/6  0  0  0  0 

2/3  1/6  0  0  0 

1/6  2/3  1/6  0  0 

0  1/6  2/3  1/6  0 

0  0  1/6  2/3  1/6 

0  0  0  1/6  2/3 


(125) 
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The  derivative  projection  matrix  D,  which  recall  is  the  same  as  the  differentiation 
matrix  V  for  periodic  boundary  conditions,  is 

D  =  B~lA.  (126) 

A  10  by  10  version  of  the  derivative  projection  matrix  for  the  current  scenario  is, 


D  = 


0 

.804 

-.215 

.057 

-.014 

0 

.014 

-.057 

.215 

-.804 

-.804 

0 

.804 

-.215 

.057 

-.014 

0 

.014 

-.057 

.215 

.215 

-.804 

0 

.804 

-.215 

.057 

-.014 

0 

.014 

-.057 

-.057 

.215 

-.804 

0 

.804 

-.215 

.057 

-.014 

0 

.014 

.014 

-.057 

.215 

-.804 

0 

.804 

-.215 

.057 

-.014 

0 

0 

'.014 

-.057 

.215 

-.804 

0 

.804 

-.215 

.057 

-.014 

-.014 

0 

.014 

-.057 

.215 

-.804 

0 

.804 

-.215 

.057 

.057 

-.014 

0 

.014 

-.057 

.215 

-.804 

0 

.804 

-.215 

-.215 

.057 

-.014 

0 

.014 

-.057 

.215 

-.804 

0 

.804 

.804 

-.215 

.057 

-.014 

0 

.014 

-.057 

.215 

-.804 

0 

(127) 

6.2  The  Differentiation  Matrix  for  a  3-rd  Order  B-Spline 

The  3-rd  order  B-spline  is  the  convolution  of  two  1-st  order  B-splines.  The  (non- 
normalized)  3-rd  order  B-spline  is  defined  as,  fl3{x)  =  (z  +  2)3  for  x  €  [-2,-1), 
/?3(x)  =  4  —  6x2  —  3x3  for  x  €  [—1,0),  /?3(x)  =  4  —  6x2  +  3x3  for  x  €  [0,1), 
/?3(x)  =  (2  —  x)3  for  x  €  [1,2],  and  0  otherwise. 


6.2.1  The  Filters  a,  6,  and  d 

The  explicitly  calculated  filters  for  the  3-rd  order  B-spline  are, 

o  =  {-1/20,  -14/5,  -49/4, 0, 49/4, 14/5, 1/20}, 

and 


(128) 


V  =  {1/140, 6/7, 1.191/140, 604/35, 1191/140,6/7, 1/140}.  (129) 

Normalization  is  not  necessary  since  the  filter  a  is  divided  by  the  filter  67.  The 
semi-discrete  Fourier  transforms  are, 


a(£)  =  *(49/2sin(£)  +  28/5sin(2£)  +  i/10sin(3£))> 


(130) 
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and 


67(0  =  604/35  +  1191/70  cos(£)  +  12/7cos(2^)  +  l/70cos(3^).  (131) 


Again,  the  semi-discrete  Fourier  transforms  are  related  by, 


WK)  =  ««)• 


(132) 


In  a  manner  analogous  to  the  previous  subsection,  we  equate  powers  in  the  Taylor 
series  of  the  semi-discrete  Fourier  transforms  and  come  to  the  conclusion  that, 


d(0  =  i£  +  cf  +  hot. 


(133) 


This  is  precisely  the  condition  that  is  required  for  8-th  order  accuracy. 


6.2.2  The  Differentiation  Matrix 


The  matrices  for  the  filters  a  and  b  are, 


A 


'  0 

49/4 

14/5 

1/20 

0 

-49/4 

0 

49/4 

14/5 

1/20 

-14/5 

-49/4 

0 

49/4 

14/5 

-1/20 

-14/5 

-49/4 

0 

49/4 

0 

-1/20 

-14/5 

-49/4 

0 

0 

0 

-1/20 

-14/5 

-49/4 

0 

0 

0 

-1/20 

-14/5 

1/20 

0 

0 

0 

-1/20 

14/5 

1/20 

0 

0 

0 

49/4 

14/5 

1/20 

0 

0 

and 


(134) 


0 

0 

-1/20 

-14/5 

-49/4  ' 

0 

0 

6 

-1/20 

-14/5 

1/20 

0 

0 

0 

-1/20 

14/5 

1/20 

0 

0 

0 

49/4 

14/5 

1/20 

0 

0 

0 

49/4 

14/5 

1/20 

0 

-49/4 

0 

49/4 

14/5 

1/20 

-14/5 

-49/4 

0 

49/4 

14/5 

-1/20 

-14/5 

-49/4 

0 

49/4 

0 

-1/20 

-14/5 

-49/4 

0 

B  =  (135) 


'  604/35 

1191/140 

6/7 

1/140 

0 

0 

1/140 

6/7 

1191/140  ' 

1191/140 

604/35 

1191/140 

6/7 

1/140 

0 

0 

1/140 

6/7 

6/7 

1191/140 

604/35 

1191/140 

6/7 

1/140 

0 

0 

1/140 

1/140 

6/7 

1191/140 

604/35 

1191/140 

6/7 

1/140 

0 

0 

0 

1/140 

6/7 

1191/140 

604/35 

1191/140 

6/7 

1/140 

0 

0 

0 

1/140 

6/7 

1191/140 

604/35 

1191/140 

6/7 

1/140 

1/140 

0 

0 

1/140 

6/7 

1191/140 

604/35 

1191/140 

6/7 

6/7 

1/140 

0 

0 

1/140 

6/7 

1191/140 

604/35 

1191/140 

.  1191/140 

6/7 

1/140 

0 

0 

1/140 

6/7 

1191/140 

604/35 
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The  derivative  projection  matrix  is 


D  =  B-'A , 

D  =  (136) 


0 

.938 

-.387 

.184 

-.078 

0 

.078 

-.184 

.387 

-.938 

-.938 

0 

.938 

-.387 

.184 

-.078 

0 

.078 

-.184 

.387 

.387 

-.938 

0 

.938 

-.387 

.184 

-.078 

0 

.078 

-.184 

-.184 

.387 

-.938 

0 

.938 

-.387 

.184 

-.078 

0 

.078 

.078 

-.184 

.387 

-.938 

0 

.938 

-.387 

.184 

-.078 

0 

0 

.078 

-.184 

.387 

-.938 

0 

.938 

-.387 

.184 

-.078 

-.078 

0 

.078 

-.184 

.387 

-.938 

0 

.938 

-.387 

.184 

.184 

-.078 

0 

.078 

-.184 

.387 

-.938 

0 

.938 

-.387 

-.387 

.184 

-.078 

0 

.078 

-.184 

.387 

-.938 

0 

.938 

.938 

-.387 

.184 

-.078 

0 

.078 

-.184 

.387 

-.938 

0 
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7  Nonperiodic  Boundary  Conditions 


In  this  section  the  differentiation  matrix  will  be  given  for  one  particular  boundary 
construction.  The  construction  considered  here  is  one  of  truncated  B-splines.  That 
is,  away  from  the  boundaries  one  simply  shifts  the  B-splines  by  one  to  generate  a 
spline  basis.  At  the  boundaries  we  simply  continue  this  shifting.  If  part  of  the  B- 
spline  goes  beyond  the  domain  then  simply  truncate.  This  construction  maintains 
the  approximation  properties  of  the  spline  subspace.  That  is,  any  n-th  order  spline 
can  be  generated  across  the  domain  using  n-th  order  B-splines  truncated  in  this  way. 


7.1  1st  Order  Spline 


Using  the  boundary  construction  outlined  in  the  previous  paragraph  and  the  1st  order 
B-spline  one  gets  the  following  matrices: 
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and  the  matrix  B  is, 
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Previously  it  was  shown  that  the  differentiation  matrix  has  the  form, 
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Table  1:  Calculated  order  of  accuracy  for  a  1st  order  spline- based  wavelet  differenti¬ 
ation  matrix  for  truncated  B-spline  boundary  construction. 

For  n  =  1  the  matrix  6  is  the  identity. 

7.1.1  Computed  Accuracy 

Again,  one  can  estimate  the  order  of  accuracy  by  applying  the  above  defined  matrix 
T>  to  evenly-spaced  samples  of  the  function 

1 

2  +  sin(x)’ 

and  comparing  to  the  exact  derivative. 

In  table  (7.1.1)  ‘Inner  Error’  denotes  the  error  at  the  middle  50%  of  the  grid 
points,  and  the  boundary  error  is  the  error  at  the  two  outermost  grid  points  at  each 
end  of  the  interval.  From  the  table  it  can  be  seen  that  the  superconvergence  is  lost  at 
the  boundaries.  That  is,  in  the  middle  of  the  interval  the  4th-order  superconvergence 
is  obtained,  whereas  at  the  boundary  the  differentiation  is  only  lst-order  accurate. 

7.2  Example  for  3-rd  Order  Spline 

For  the  3rd  order  B-splines  the  calculations  are  similarly  straightforward.  There  is, 
however,  one  problem:  the  matrix  6  is  ill-conditioned.  We  need  the  inverse  of  b  and 
for  grids  larger  than  about  32  points  the  errors  introduced  by  inverting  b  make  it 
very  difficult  to  determine  the  accuracy.  One  can  avoid  this  problem  when  testing 
accuracy  by  simply  not  including  b  in  the  calculations.  That  is,  b  comes  from  the 
quadrature  equations  and  plays  no  role  in  the  differencing.  The  following  are  the 
matrices  A  and  B  produced  by  the  truncated  B-spline  boundary  construction: 
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7.2.1  Computed  Accuracy 

If  we  apply  the  matrix  V  =  B~l  A  to  the  evenly-spaced  samples  of  the  function 


^ X ^  2  +  sin(£x)  ’ 

for  a  fixed  wave  number  £  it  is  difficult  to  determine  the  accuracy  both  at  the  bound¬ 
aries  and  in  the  middle  of  the  interval  at  the  same  time.  That  is,  it  was  proven  that 
the  accuracy  ‘in  the  middle’  is  2n  +  2  or  8  when  n  =  3.  For  the  current  boundary 
construction,  however,  the  accuracy  is  3  at  the  boundaries.  This  vast  difference  in 
accuracy  requires  that  the  frequency  of  the  test  function  be  low  in  order  to  test  the 
boundary  accuracy  and  high  in  order  to  test  the  superconvergence  in  the  middle  of 
the  interval. 

In  table  (7.2.1)  the  wave  number  is  £  =  2x.  This  wave  number  is  sufficiently  high 
to  test  the  superconvergence. 


V  —  B~XA 
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Table  2:  Calculated  order  of  accuracy  for  a  3rd  order  spline-based  wavelet  differen¬ 
tiation  matrix  for  the  truncated  B-spline  boundary  construction.  The  test  function 
has  a  high  frequency  in  order  to  capture  the  8th  order  differentiation  accuracy. 
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Table  3:  Calculated  order  of  accuracy  for  a  3rd  order  spline-based  wavelet  differen¬ 
tiation  matrix  for  the  truncated  B-spline  boundary  construction.  The  test  function 
has  a  low  frequency  in  order  to  capture  the  3rd  order  differentiation  accuracy. 

Now  the  lower  frequency  £  =  1  is  used  to  test  the  accuracy  at  the  boundaries.  See 
table  (7.2.1).  Note  that  in  table  (7.2.1)  that  a  few  of  the  entries  in  the  ‘Error  Ratio’ 
columns  are  blank.  This  is  done  whenever  the  ‘Inner  Error’  is  of  the  order  10~1S  and 
can  no  longer  be  resolved  on  the  computer.  Table  (7.2.1)  indicates  the  accuracy  only 
at  the  boundaries.  As  in  table  (7.2.1)  it  can  be  seen  that  the  boundary  accuracy  is 
3rd  order  and,  consequently,  the  superconvergence  is  lost  at  the  boundaries. 


The  next  section  contains  a  few  comments  on  this  loss  of  superconvergence. 


8  Loss  of  Superconvergence 


In  [5]  it  was  shown  that  under  a  given  boundary  construction  the  Daubechies- based 
wavelet  superconvergence  which  was  proved  in  [4]  is  lost.  In  this  paper  we  have 
encountered  a  similar  phenomenon.  I  believe  that  this  in  not  simply  a  function  of 
the  the  particular  boundary  functions  which  were  chosen  but  a  general  characteristic 
for  numerical  methods  which  are  based  on  approximation  by  translation.  In  fact, 
Gottlieb,  et.  al,  see  [3],  have  shown  that  the  superconvergence  encountered  in  the 
finite  element  method  under  periodic  boundary  conditions  cannot  be  recovered  when 
nonperiodic  boundary  conditions  are  present:  '...  no  matter  how  the  approximating 
subspace  Sk  is  chosen,  the  superconvergence  property  is  lost  if  there  are  characteristics 
leaving  the  domain.’ 

Superconvergence  is  a  property  which  depends  on  the  overlap  of  bases  functions 
with  neighboring  bases  functions.  At  the  boundaries  the  overlap  only  occurs  on  one 
side.  For  a  Daubechies- based  wavelet  system  where  the  bases  functions  are  orthogonal 
there  is  a  minimum  amount  of  overlap  which  must  occur  at  the  boundary  in  order 
to  maintain  the  superconvergence.  That  is,  for  a  Daubechies  4  wavelet,  4-th  order 
differentiation  accuracy  is  encountered.  In  order  to  maintain  the  superconvergence 
one  must  have  at  least  3rd  order  differentiation  accuracy  at  the  boundaries.  That  is, 
four  conditions  must  be  satisfied  requiring  an  overlap  of  four  basis  functions  at  the 
boundary.  As  of  now,  there  are  no  boundary  constructions  satisfying  this  criterion. 

For  spline-based  wavelet  systems  the  situation  is  similar  but  slightly  more  compli¬ 
cated  since  the  underlying  functions,  B-splines,  are  not  orthogonal.  Again,  however, 
a  general  rule-of-thumb  is  that  the  overlap  of  the  B-splines  at  the  boundary  should 
be  the  same  as  the  superconvergence  encountered  with  periodic  boundary  conditions. 
That  is,  for  the  1st  order  B-spline  we  have  shown  here  that  the  differentiation  is  4th 
order  accurate  where  periodic  boundary  conditions  are  assumed.  This  would  require 
an  overlap  of  4  boundary  functions  in  order  to  have  3rd  order  boundary  accuracy 
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which  maintains  the  4th  order  overall  accuracy. 
The  next  section  concludes  this  paper. 
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9  Conclusion 


Let  me  restate  a  few  of  the  salient  features  of  this  paper.  First  of  all,  circularity  plays 
a  prominent  role  when  periodic  boundary  conditions  are  assumed.  The  benefit  of 
working  with  circular  matrices  is  that  they  commute.  In  the  periodic  case  this  implies 
that  the  differentiation  matrix  is  the  same  as  the  derivative  projection  matrix, 

V  =  C~lDC  =  D. 

The  proof  of  superconvergence  is,  however,  the  most  meaningful  result  of  the  pa¬ 
per.  I  have  shown  that  the  matrix  D  displays  accuracy  of  order  2n  +  2  where  n  is 
the  order  of  the  underlying  B-spline  used  to  construct  the  scaling  functions.  Also,  it 
was  seen  that  this  accuracy  holds  regardless  of  of  whether  the  scaling  functions  are 
chosen  to  be  orthogonal  under  translation  or  not.  That  is,  the  scaling  functions  can 
be  required  to  be  B-splines,  C-splines,  or  orthogonal.  The  parameters  which  dictate 
the  properties  of  the  scaling  function  ‘divide  out’  when  the  differentiation  matrix  is 
calculated.  Because  of  this  feature,  only  one  theory  is  necessary  to  illustrate  the  ac¬ 
curacy  for  the  spline-based  differentiation  matrix  when  periodic  boundary  conditions 
are  assumed. 

When  the  boundary  conditions  are,  on  the  other  hand,  not  periodic  it  has  been 
shown  that  for  boundary  functions  constructed  from  truncated  B-splines  that  the  su¬ 
perconvergence  is  lost  at  the  boundary.  I  speculate  that  this  loss  of  superconvergence 
will  continue  to  be  a  characteristic  of  wavelet  bases  defined  on  an  interval. 

It  is,  also,  interesting  to  note  the  relationship  between  spline-based  wavelet  meth¬ 
ods  and  compact  finite  difference  methods.  That  is,  the  underlying  character  of  a 
spline-based  wavelet  method  can  be  understood  in  terms  of  compact  finite  difference 
schemes. 
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